UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADP014810 

TITLE:  A Self-Contained  Mapping  Closure  Approximation  for  Scalar 
Mixing 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  Annual  Research  Briefs  - 2003  [Center  for  Turbulence  Research] 
To  order  the  complete  compilation  report,  use:  ADA420749 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADPO 14788  thru  ADPO 14827 


UNCLASSIFIED 


Center  for  Turbulence  Research 
Annual  Research  Briefs  2003 


277 


A self-contained  mapping  closure  approximation 

for  scalar  mixing 

By  Guo- Wei  He  f and  Zi-Fan  Zhang  | 


1.  Objective  and  motivation 

Scalar  turbulence  exhibits  interplays  of  coherent  structures  and  random  fluctuations 
over  a broad  range  of  spatial  and  temporal  scales.  This  feature  necessitates  a probabilistic 
description  of  the  scalar  dynamics,  which  can  be  achieved  comprehensively  by  using  prob- 
ability density  functions  (PDFs).  Therefore,  the  challenge  is  to  obtain  the  scalar  PDFs 
(Lundgren  1967;  Dopazo  1979).  Generally,  the  evolution  of  a scalar  is  governed  by  three 
dynamical  processes:  advection,  diffusion  and  reaction.  In  a PDF  approach  (Pope  1985), 
the  advection  and  reaction  can  be  treated  exactly  but  the  effect  of  molecular  diffusion 
has  to  be  modeled.  It  has  been  shown  (Pope  1985)  that  the  effect  of  molecular  diffusion 
can  be  expressed  as  conditional  dissipation  rates  or  conditional  diffusions.  The  currently 
used  models  for  the  conditional  dissipation  rates  and  conditional  diffusions  (Pope  1991) 
have  resisted  deduction  from  the  fundamental  equations  and  are  unable  to  yield  satisfac- 
tory results  for  the  basic  test  cases  of  decaying  scalars  in  isotropic  turbulence,  although 
they  have  achieved  some  success  in  a variety  of  individual  cases.  The  recently  developed 
mapping  closure  approach  (Pope  1991;  Chen,  Chen  & Kraichnan  1989;  Kraichnan  1990; 
Klimenko  & Pope  2003)  provides  a deductive  method  for  conditional  dissipation  rates 
and  conditional  diffusions,  and  the  models  obtained  can  successfully  describe  the  shape 
relaxation  of  the  scalar  PDF  from  an  initial  double  delta  distribution  to  a Gaussian  one. 
However,  the  mapping  closure  approach  is  not  able  to  provide  the  rate  at  which  the  scalar 
evolves.  The  evolution  rate  has  to  be  modeled.  Therefore,  the  mapping  closure  approach 
is  not  closed.  In  this  Letter,  we  will  address  this  problem. 

The  evolution  rate  of  scalar  is  a key  quantity  in  modeling  turbulent  mixing  for  both 
conserve  and  reactive  scalars  (Cha  & Trouillet  2003).  It  specifies  the  characteristic  time 
scale  of  scalar  evolution.  It  has  been  shown  that  the  decay  rate  of  scalar  depends  on  the 
relative  length  scale  ratio  of  the  initial  scalar  and  velocity  fields  (Warhaft  & Lumley  1978; 
Sreenivasan,  Tavoularis,  Henry  & Corrsin  1980;  Durbin  1982;  Eswaran  & Pope  1988; 
Mell,  Kosaly  & Riley  1991),  and  recently,  the  asymptotic  decay  of  scalar  turbulence 
has  been  extensively  studied  (Eyink  & Xin  2000;  Chaves,  Eyink,  Frisch  & Vergassola 
2001;  Chertkov  & Lebedev  2003).  Nearly  all  existing  models  for  scalar  mixing,  ranging 
from  the  simple  (conditional)  moment  approaches  to  the  full  PDF  approaches,  require 
information  on  the  time  scales.  These  models  are  mainly  based  on  the  assumption  of 
a direct  proportionality  between  the  scalar  time  scales  and  the  turbulence  time  scales 
(Fox  1995).  Moreover,  they  exclude  the  effects  of  chemical  reaction  on  the  time  scales 
of  scalar  evolution.  The  mapping  closure  approach  of  time-dependent  reference  fields 
(Girimaji  1992)  can  provide  the  time  scale  externally,  which  highlights  an  attack  line  to 
this  problem. 
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We  develop  here  a mapping  closure  approximation  (MCA)  approach  for  the  time  scale 
of  scalar  evolution.  The  MCA  approach  is  based  on  multi-point  joint  PDFs.  In  the  map- 
ping closure  approach  (Pope  1991;  Chen,  Chen  & Kraichnan  1989),  an  unknown  random 
field  is  mapped  from  a known  random  Gaussian  field  so  that  the  evolution  of  the  un- 
known random  field  can  be  described  by  the  mapping  function.  The  mapping  function 
is  obtained  from  the  transport  equation  for  the  unknown  random  field  and  the  Gaussian 
closure.  Since  the  mapping  function  is  constructed  at  the  level  of  one-point  PDFs,  it  is 
not  able  to  provide  the  information  on  two-point  statistics,  such  as  the  time  scales  (He, 
Rubinstein  & Wang  2002).  In  the  MCA  approach,  the  mapping  functions  are  constructed 
at  the  levels  of  multi-point  joint  PDFs.  The  mapping  function  based  on  the  two-point 
joint  PDFs  could  provide  the  correct  information  on  time  scales. 


2.  Main  results 

We  consider  the  simple  case  of  a reactive  scalar  advected  by  a stochastic  velocity  field: 

^ + u-v¥3  = r vV  + QM,  (2.i) 

where  the  velocity  field  u obeys  V ■ u = 0 and  is  independent  of  the  scalar  field.  With- 
out the  loss  of  generality,  it  may  be  prescribed  as  a known  homogeneous  and  isotropic 
Gaussian  field.  T is  a molecular  diffusivity,  Q{<p)  mimics  a one-species  chemical  reaction. 

In  the  MCA  approach,  the  scalar  field  is  mapped  from  a known  random  field  by  a 
mapping  function 


<p(x,t)  = X(0(x,t),t).  (2.2) 

Here,  the  known  random  field  6(x,  t ) is  taken  as  a Gaussian  reference  field.  Its  one-point 
and  two-point  joint  PDFs  are  defined  by 


21 


9i(v)  = -7==  exp 

V27T 

92{Vi,V2,r,t)  =g(m,V2,p(r,t))  = 


exp 


where 


p(r,  t)  = 


2-Ky/l  — p2 
(6(x,t)6(x  + r,t) 


Vi  + V2  ~ 2PV1V2 


2(1  - p2) 


(2.3) 
, (2.4) 


(2.5) 


<02(x,t)) 

Here  r is  the  magnitude  of  separation  vector  r.  The  mapping  function  is  required  to  rep- 
resent the  one-point  and  two-point  joint  PDFs  of  the  scalar  <p  via  the  following  equations 


=9i{v) 


dX(V,t) 


drj 


-l 


f2(ipi,tp2,r,t)  = 92{VuV2,r,t) 


(2.6) 

(2.7) 


-dX(m,t)  dXjrj^t)]-1 

dr)i  dr)2 

In  the  classic  mapping  closure  approach  (Pope  1991;  Chen,  Chen  & Kraichnan  1989), 
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the  mapping  function  X is  only  required  to  represent  the  one-point  PDF  of  the  scalar  <p 
via  equation  (2.6).  Differentiating  equation  (2.6)  with  respect  to  t yields 


dfi  d_  r dx' 

dt  dip  1dt 


= 0. 


(2.8) 


Meanwhile,  the  transport  equation  for  the  one-point  PDF,  fi(ip,t),  can  be  derived  by 
the  test  function  method  (Gotoh  & Kraichnan  1993;  Kimura  & Kraichnan  1993)  as 


a-k  + ^^  + QW)\^  ]=°-  (2-9) 

Thus,  comparing  the  two  equations  (2.8)  and  (2.9)  with  the  substitution  of  (2.2),  we 
obtain  the  exact  result 


dX 

dt 


= r(Wl  <P  = X(r,,t))  + Q(<p). 


(2.10) 


The  conditional  moment  in  (2.10)  can  be  evaluated  from  the  mapping  function  (2.2)  and 
the  Gaussianity  (2.3)  and  (2.4)  of  the  reference  field  9 (Panchev  1971) 


r(v2<p\<P  = x(T1,t))  = -Cp"(o,t)r 


o2x  ax 


dr]2  ^ dp  J ’ 

where  the  prime  denotes  the  derivative  with  respect  to  separation  r.  As  a result,  the 
transport  equation  for  the  mapping  function  (2.2)  becomes 


(2.11) 


dX 

dt 


-Cp"(0,t)T 


d2x  ax 


[ dr]2  71  dr]  J 

It  is  easily  shown  from  the  Gaussianity  (2.3)  and  (2.4)  that  ((V#)2)  = — Cp"(0,t ), 
where  C = 2 for  a two-dimensional  physical  space  and  C — 3 for  a three-dimensional 
physical  space.  Equation  (2.12)  has  been  obtained  in  (Chen,  Chen  & Kraichnan  1989), 
where  — Cp  (0,  t)  is  represented  by  the  variance  ((W)2).  However,  the  correlation  p(r,  t) 
in  equation  (2.12)  still  remains  to  be  unknown  and  has  to  be  input  externally.  For  ex- 
ample, it  is  set  using  the  results  from  direct  numerical  simulation  in  (Chen,  Chen  & 
Kraichnan  1989).  Therefore,  equation  (2.12)  is  unclosed. 

The  two-point  correlation  p(r,t)  cannot  be  obtained  from  the  one-point  PDF  gi(r]). 
Rather,  it  has  to  be  calculated  from  the  two-point  joint  PDF  92(91,92,  r,t).  Hence,  we 
propose  to  invoke  the  two-point  joint  PDF  (2.4),  which  is  not  used  in  the  classic  mapping 
closure  approach.  By  differentiating  (2.7)  with  respect  to  t,  we  obtain 


+ Q(X). 


(2.12) 


0/a  + d 


h 


dXx 


dt  dip  1 [ dt  J dip2  1 dt 


+ 


d 


/2 


dX 2 


h dg2 

92  dt  ' 


(2.13) 


The  transport  equation  for  the  two-point  joint  PDF  fz  derived  from  the  test  function 
method  (Gotoh  & Kraichnan  1993;  Kimura  & Kraichnan  1993)  has  the  form 


^£  + Vr-[/2((u2-U1)|V>1^2)]  = 

-^[/2(rvV +<5(91)1^1^2)]  -^[/2(rvV2  + ^2)|V'i,^2)].  (2.14) 
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Subtracting  (2.13)  from  (2.14)  leads  to 

+ vr  ■ [/2((U2  - Ul)  I V’l.V’a)]  = (2.15) 

where 

Hk  = r<vV*  I A)  - r (vV*  I Vb,  V>2).  (2.16) 

The  first  term  in  (2.16)  is  the  conditional  diffusion  on  a single  given  scalar  and  is 
calculated  in  equation  (2.11).  The  second  term  is  the  conditional  diffusion  on  two  given 
scalars  at  two  different  locations  and,  again,  can  be  evaluated  using  the  mapping  function 
(2.2)  and  the  Gaussianity  (2.3)  and  (2.4)  of  the  reference  field  9 (Panchev  1971) 


(V2<Pk\A,A)  = ^^<(V6>fc)2jj?i,7/2)  + < V26>fc|  771,772) 


= {~cp"M + (i-Mlw + - p{r^2  - 1 + p2(r^ } 


+ 


dXk 


dr)k  1 -p2{r,t)  [ 


p"(r,  t)  + - ) (vi  ~ P(r,t)r}k ) + Cp"(0,t)(r)k  - p{r,t)rn) 


(?.17) 


where  l = 1 for  k — 2 and  l = 2 for  k = 1.  Multiplying  (2.15)  by  ip  1 and  ip2  and 
then  taking  the  mean  with  substitution  of  (2.2),  (2.4),  (2.7)  and  (2.17),  we  obtain  the 
transport  equation  for  p(r,  t)  as  follows 


p"(r,  t)  + - - Cp(r,  t)p"( 0,  t)  + p,2(r,  t ) 


-1 


= 2r- 
d2Xx  d2X2 


dm  dm  / \ dm  dm 


dXx  dX2 


2.18) 


Equations  (2.12)  and  (2.18)  form  a closed  system  for  the  mapping  function,  where 
equation  (2.12)  describes  the  evolution  of  the  shape  of  the  mapping  function  and  (2.18) 
specifies  the  rate  at  which  the  mapping  function  evolves.  In  equation  (2.18),  the  second 
term  on  the  left-hand  side  corresponds  to  advection,  the  first  three  terms  on  the  right- 
hand  side  correspond  to  diffusion  and  the  last  term  on  the  right-hand  side  corresponds  to 
the  effect  of  nonlinear  mapping.  The  last  term  vanishes  if  the  mapping  function  is  linear. 
We  note  that  p(r,  t)  is  the  correlation  function  of  the  reference  field  and  is  dependent  on 
the  mapping  function. 

The  realizability  condition  of  equation  (2.18)  is  \p{r,t)\  < 1.  For  the  pure  diffusion 
processes  with  the  initial  Gaussian  distributions  of  positive  correlations,  the  diffusion 
term  p (r,  t)  + p (r,  t)/r  and  the  damping  term  -Cp"  (0,  t)p(r,  t)  decrease  the  amplitudes 
of  the  correlation  p(r,t),  so  that  the  solution  of  equation  (2.18)  is  realizable. 

The  mapping  equation  (2.12)  is  closed  using  the  two- point  joint  PDF  constraint  (2.7), 
from  which  the  correlation  equation  (2.18)  is  derived.  Another  possibility  for  the  clo- 
sure is  to  use  the  constraint  of  the  joint  PDF  for  the  scalar  and  its  derivative  (Chen, 
Chen  &;  Kraichnan  1989),  which  leads  to  an  unclosed  equation  for  p(0,  t)  and  its  spatial 
derivatives.  It  points  to  another  direction  to  go  beyond  the  one-point  mapping  (2.2)  for 
different  purpose. 

The  performance  of  the  MCA  models  (2.12)  and  (2.18)  are  evaluated  against  direct  nu- 
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Figure  1.  Pure  diffusion  for  the  conserved  scalar  with  a double-delta  initial  field.  The  left: 
variance  and  dissipation  vs  time.  Solid  lines  for  MCA  and  squares  for  DNS;  The  right:  probability 
density  function.  Solid,  dashed  and  dash-dotted  lines  for  MCA  at  time  0.05,0.15  and  0.5,  and 
circles,  squares  and  deltas  for  DNS  at  the  same  time,  respectively. 


merical  simulation.  The  scalar  equation  (2.1)  and  the  MCA  modeled  system  of  equations 
(2.12)  and  (2.18)  are  numerically  solved  in  a cyclic  square  of  side  27r,  using  second- 
order  Adams-Bashforth  scheme  in  time  and  fourth-order  central  finite-difference  scheme 
in  space.  In  all  the  cases,  the  non-dimensional  molecular  diffusivity  T = 0.01.  Boundary 
conditions  are  periodic  in  space,  except  that  the  ones  in  the  direction  of  the  reference  field 
are  obtained  by  extrapolation.  The  initial  fields  for  (2.1)  are  double-delta  distributions 
or  isotropic  Gaussian  distributions  with  their  energy  spectra  Ev(k)  oc  fc~17/3.  Thus,  the 
initial  mapping  for  (2.12)  and  the  initial  correlation  for  (2.18)  can  be  calculated  from 
their  definitions. 

The  velocity  field  is  a given  homogeneous  isotropic  Gaussian  process  with  spectrum  of 
the  form  E(k)  oc  fc~5/3.  We  will  here  present  the  results  for  the  rapidly  changed  velocity. 
That  is,  the  velocity  fields  are  set  to  change  at  each  time  step  ( Kimura  & Kraichnan 
1993).  In  order  to  isolate  the  effects  of  the  MCA  models  on  diffusivity  and  reaction,  the 
advection  term  in  equation  (2.18)  are  calculated  directly  from  the  DNS  without  invoking 
any  models.  The  analytical  treatment  on  this  term  can  be  found  in  ( Kimura  & Kraichnan 
1993). 

Figure  1 compares  the  evolutions  of  the  variance  (ip2(x,t))  and  the  PDF  fi{ip,t)  ob- 
tained by  DNS  of  equation  (2.1)  with  an  integration  of  the  MCA  models  (2.12)  and  (2.18) 
for  the  case  of  pure  diffusion:  u = 0 and  Q = 0.  The  initial  condition  is  set  as  a double 
delta  distribution.  The  results  illustrate  that  the  MCA  models  represent  the  relaxation  of 
the  double  delta  PDF  to  the  Gaussian  PDF  not  only  in  its  shape  but  also  at  the  correct 
rate  of  evolutions. 

In  Fig.  2,  the  same  plots  are  made  for  the  diffusion-reaction  equation:  u = 0 and 
Q(<p)  = — 20<p\(p\,  with  the  initial  Gaussian  distribution.  The  MCA  models  are  also  in 
agreement  with  the  DNS  results.  It  shows  that  the  MCA  models  can  represent  the  effects 
of  both  diffusion  and  reactions. 

Further  comparisons  are  made  in  Fig.  3 for  the  advection-diffusion  equation  with  the 
same  initial  condition  as  used  in  Fig.  2,  where  the  Peclet  number  is  about  101.  Evidently 
the  MCA  models  are  in  good  agreement  with  the  decay  rate  of  the  scalar. 


Figure  3.  Advection-diffusion  for  the  initial  Gaussian  scalar.  The  left:  variance  and  dissipation 
vs  time.  Solid  lines  for  MCA  and  squares  for  DNS;  The  right:  probability  density  function.  Solid, 
dashed  and  dash-dotted  lines  for  MCA  at  time  0.0, 0.05  and  0.3,  and  circles,  squares  and  deltas 
for  DNS  at  the  same  time,  respectively. 


In  this  Letter,  the  two-point  joint  PDF  is  represented  by  a known  two-point  joint 
Gaussian  PDF  and  a mapping  function.  The  correlation  of  the  joint  Gaussian  PDF  and 
the  mapping  function  evolve  under  their  dynamics.  Therefore,  the  representation  of  the 
two-point  joint  PDF  evolves  with  the  dynamics  of  equation  (2.1)  via  the  correlation 
and  mapping  function.  If  it  happens  that  p{r,t)  = 0,  equation  (2.7)  is  simplified  to 

This  is  the  stochastic  ansatz  in  the  BBGKY  hierarchy  in 
statistical  mechanics  (Balescu  1975),  where  the  correlations  are  specified  a priori  and  then 
fixed.  The  MCA  approach  does  not  invoke  this  kind  of  a priori  assumption.  Moreover, 
a general  IV-point  (N  > 2)  joint  PDF  can  also  be  obtained  using  the  known  two-point 
joint  Gaussian  PDF  and  the  mapping  function  (2.2),  since  it  can  be  similarly  represented 


283 


A self-contained  mapping  closure  approximation 

by  the  mapping  function  (2.2)  and  the  iV-point  joint  Gaussian  PDF.  The  latter  is  fully 
determined  by  its  two-point  joint  PDF.  In  this  sense,  the  MCA  approach  is  self-contained. 


3.  Conclusion  and  future  work 

A self-contained  MCA  approach  is  developed  for  modeling  scalar  mixing  in  a stochas- 
tic velocity  field.  It  differs  from  the  classic  mapping  closure  approach  in  that  the  MCA 
approach  makes  use  of  two-point  PDFs  to  represent  the  time-evolving  correlations  of  the 
reference  fields  and  thus  the  scalar  fields.  Unlike  usual  treatments  in  the  BBGKY  hierar- 
chy (Balescu  1975),  where  the  representations  are  specified  a priori,  the  representations  in 
the  MCA  approach  are  allowed  to  evolve  in  coordinate  with  the  dynamics  of  scalar  mixing. 
The  results  obtained  using  the  new  approach  are  in  good  agreement  with  the  DNS  results 
for  the  three  cases  of  pure  diffusion,  diffusion-reaction  and  advection-diffusion.  The  ap- 
proach is  under  further  development  for  more  complex  situations  including  multi-scalar 
mixing  and  inhomogeneous  scalar  fields,  using  time-evolving  Gaussian  or  non-Gaussian 
reference  fields. 
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